Mechanism of High Temperature Superconductivity in a striped Hubbard Model 
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It is shown, using asymptotically exact methods, that the two dimensional repulsive Hubbard 
model with strongly modulated interactions exhibits "high temperature superconductivity" . Specif- 
ically, the explicit modulation, which has the same symmetry as period 4 bond-centered stripes, 
breaks the system into an alternating array of more and less heavily hole doped, nearly decoupled 
two-leg ladders. It is shown that this system exhibits a pairing scale determined by the spin-gap of 
the undoped two-leg ladder, and a phase ordering temperature proportional to a low positive power 
of the inter-ladder coupling. 
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Much has been written concerning the mechanism of 
high temperature superconductivity (HTC) since the dis- 
cover^ of the cuprate superconductors in 1986, and in- 
deed even before that. However, what is meant by "the 
mechanism" is rarely defined, and clearly evokes differ- 
ent images for different authors. The BCS mechanism, 
in which pairing is a consequence of a weak induced at- 
traction produced by the exchange of phonons between 
well defined quasiparticles, is not only consistent with a 
remarkable number of experimental facts in conventional 
superconductors, it is also of well established theoretical 
validity in simple models. Because it is a weak coupling 
theory, even the mean-field estimate of T c (which is expo- 
nentially small, T c oc exp(— 1/g) where g is the induced 
attraction) is known to be quantitatively reliable 2 . How- 
ever, there are many well known reasons to believe that 
the BCS mechanism always leads to low T c 's as recently 
reviewed in Ref. 0. 

An alternative idea, which has been the focus of 
much of the theoretical effort in the field, is that in a 
doped Mott insulator, high temperature superconductiv- 
ity arises directly from the repulsive interactions between 
electrons. However, even as a point of principle, the va- 
lidity of a mechanism of this sort has not been well es- 
tablished for any simple model. 

In this note we demonstrate the existence of a "high 
temperature superconducting" phase of the Hubbard and 
t — J models on a square lattice with periodically modu- 
lated parameters^ - see Eq. (0. In particular, we show 
that a period 2 modulation can produce superconductiv- 
ity with a relatively low T c in a restricted doping range, 
while a period 4 modulation produces higher critical tem- 
peratures on a broader range of doping. Specifically, we 
consider a caricature of a stripe ordered state consisting 
of a quasi-one dimensional array of two leg Hubbard lad- 
ders weakly coupled to each other with a hopping ma- 
trix element St. For a range of electron densities per 
site, < n >= 1 — x, it has been well establishe d 5 ' 6 ' 7 i 8 i 9 i 10 
that the two leg ladder exhibits a Luther- Emery liquid 11 
phase, with a large spin-gap, A s ~ J/2, and a divergent 



superconducting susceptibility for T <C A s 



Xsc(T) ~ A s / T 



2-R- 1 



(1) 



where K is the charge Luttinger parameter, and T is the 
temperature. This sounds like a promising start. How- 
ever, a non-zero T c is impossible in one dimension (ID), 
so to have a chance of a high transition temperature, 
inter-ladder couplings must be taken into account. If 
all the ladders are equivalent (a caricature of a period 
2 stripe ordered or column stat o 12 i 13 ), we shall see that 
this coupling leads to a superconducting state in a re- 
stricted range of small x with rather low T c . For more 
substantial values of x, it inevitably leads to an insulat- 
ing, incommensurate charge density wave (CDW) state 
with (in units in which the lattice constant is a = 1) 
an ordering wave number P = 2ttx. (It is customary 
to call this the 4fcp-CDW since, despite the fact that 
there is a spin-gap and hence no Fermi surface whatso- 
ever, P = ifcFt where fcp is the Fermi momentum of a 
one-dimensional non-interacting electron gas at the same 
electron density.) That the superconducting transition is 
so easily preempted by CDW order follows from the fact 
that the CDW susceptibility of the Luther-Emery liquid 
diverges as 



Xcdw(P,T) ~ A, I T 



2-K 



(2) 



Under most circumstances for repulsive interactions K < 
1, and hence xcdw of Eq. J2J 1S more strongly divergent 
than xsc of Eq. ijljl. However, if we consider an al- 
ternating array of A and B type ladders (with different 
electron affinities) then the tendency to CDW order is 
greatly suppressed due to the mismatch between order- 
ing vectors, Pa and Pb, on neighboring ladders^*i^. We 
shall show that, so long as the exponent inequalities 



2 > K 



K- 



K A ; 2 > K- 



(3) 



are satisfied, the superconducting instability wins out. 
(If the Luttinger parameter is the same for both ladders, 
these inequalities reduce to K > K c = (yo — 1) w 0.8.) 
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Under these circumstances, the superconducting 
(Kosterlitz-Thouless) transition temperature can be re- 
liably estimated by treating the ID fluctuations exactly 
but the inter-ladder Josephson coupling J in mean- field 
approximationi&ii. 

A ( L\° 1K A K B , , 

\WJ 1 a ~ [4K A K B -K A -K B ] W 

where J is an effective coupling and W is a microscopic 
energy which we will discuss in detail below; typically 
we find J ~ St 2 /J and W ~ J. Although T c is small 
for small J ^ it is only power law small. In fact typi- 
cally a ~ 1. Because of the mean- field character of this 
estimate for T c , one expects this to be an upper bound 
to the actual T c . One also generally expects T c to be 
somewhat suppressed by phase fluctuations but typically 
by no more than a factor of 2. Indeed, a perturbative 
renormalization-group treatment for small J yields the 
same power law dependence as Eq.gJ suggesting that this 
expression is asymptotically exact for J << W. This 
fact is supported in Appendix [S] where the accuracy of 
interchain mean-field estimates is discussed for related 
models. 

Since we expect T c to be smooth function of St/ J , it is 
reasonable to extrapolate Eq. Q to the case in which St 
is a substantial fraction of J . This suggests a maximum 
T c of order A s , and so can easily account for relatively 
high transition temperature a 18 ' 19 . This is in contrast to 
the case of an exponentially small T c as obtained, for 
example, in a BCS-like mechanism. 

I. THE STRIPED HUBBARD MODEL 

While the results obtained in this paper are quite ro- 
bust in the sense that they apply for a broad range of 
microscopic interactions, to establish their validity it is 
useful to consider an explicit model. The model we study 
is the striped Hubbard model: 

H= - ^2 V : r'[ct<r c ?V +h.c] (5) 

<r,f /> ,(T 

where < f , r 1 > designates nearest-neighbor sites, ct 
creates an electron on site f with spin polarization 
= ±1 and satisfies canonical anticommutation rela- 
tions, and U > is the repulsion between two elec- 
trons on the same site. In the limit of strong repulsions, 
U 3> tp p , this model reduces approximately to the cor- 
responding t — J model, which operates in the subspace 
of no double occupied sites, but with an exchange cou- 
pling, Jpj« — A\tf } fi\ 2 /U between neighboring spins. Our 
results only depend on the low-energy physics of the lad- 
der and, thus, apply equally to the t — J and Hubbard 
models. 



In the translationally invariant Hubbard model, tj> t f — 
t and eft = 0. The striped version of this model is still 
translationally invariant along the stripe direction (which 
we take to be the y axis), so tp t p + y = t. However, per- 
pendicular to the stripes the hopping matrix takes on 
alternately large and small values: tpp +x = t' for r x = 
even, and = St <C t' ~ t for r x = odd. This defines 

a "period 2 striped Hubbard model," as shown in Fig. 1. 
For the "period 4 striped Hubbard model," we include a 
modulated site energy, ep = \/2e cos[7rr x /2 — ir/4], which 
has site energy e and — e respectively on every other 2-leg 
ladder, with e 3> St. 
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FIG. 1: Schematic representation of the striped Hubbard 
model analyzed in this paper. 



II. ISOLATED 2-LEG LADDER 

For St = 0, the model breaks up into a series of dis- 
connected 2-leg ladders. Considerable analytic and nu- 
merical effort has gone into studying the properties of 
2-leg t — J and Hubbard ladders, and much is known 
about them. For x = 0, the undoped two leg ladder 
has a unique, fully gapped state, referred to as C0S0 in 
the notation of Ref. 0, meaning gapless charge and 
gapless spin modes. In the large U limit, the magnitude 
of the spin-gap of the undoped-*^ ladder is A s w J/2. 
Then, for a substantial range of x (0 < x < x c ) the 
ladder exhibits a Luther-Emery or C1S0 phase, with a 
spin-gap that drops smoothlySi with increasing x, and 
vanishes at a critical value of the doping, x — x c . (This 
particular Luther-Emery liquid is knowr&Sii&il to have 
"d-wave-like" superconducting correlations, in the sense 
that the pair-field operator has opposite signs along the 
edge of the ladder (y direction) and on the rungs (x- 
direction).) For x > x c , the numerical results are scarce, 
nor is there uniform agreement concerning the number 
of phases; there may 9 - or may nofeifi be narrow ranges of 
C2S1 and C2S2 phases for x slightly larger than x c . At 
any rate, for x large enough, x c < x' c < x < 1, the lad- 
der manifestly enters a Luttinger liquid C1S1 phase, and 
finally, a trivial C0S0 phase when x = 1 (< n >= 0). 

For the purposes of the present paper, we will confine 
ourselves to the range of parameters where both A and 
B type ladders are in the Luther-Emery phase. The low 
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energy physics (at all energies less than A s ) of the two- 
leg ladder in the Luther-Emery phase is contained in the 
free bosonic Hamiltonian for the collective charge degrees 
of freedom, 



1 



K{d v ey + - (a, 



(6) 



where is the CDW phase and 9 is the supercon- 
ducting phase; these two fields are dual to each other, 
and so satisfy the canonical commutation relations, 
\<j>{y'),d y 6{y)} = i5(y — y'). This effective Hamiltonian 
is general and physical; the precise x dependence of the 
spin-gap, A s , the charge Luttinger exponent, K, the 
charge velocity, v c , and the chemical potential, fi(x), de- 
pends on details such as the values oiU/t and t'/t. For 
certain cases£*&£ these have been accurately computed in 
Monte-Carlo studies, and these studies could be straight- 
forwardly extended to other values of the parameters. 

The ellipsis in Eq. J^J represent cosine potentials, 
which we will not explicitly exhibit here, that produce 
the Mott gap A« at x = 0. A consequence of these 
terms is that for x — > 0, the elementary excitations are 
charge 2e solitons that can either be viewed as spinless 
Fermions or hard-core bosons, with a dispersion relation 
E(k) — Am + tk 2 . One consequence of this is thals^* 2 ^ 
K — * 2 and v c — > 2irtx as x — ► 0. A second consequence 
is that the renormalized harmonic theory, which retains 
only the explicitly exhibited terms in Eq. JfjJ), is valid in 
a range of energies which is small in proportion to the 
effective Fermi energy, Ep = 2-kIx 2 . (An estimate of 
t w t/2 can be obtained from the DMRG study of the 
t - J ladder with J/t = 1/3 in Ref. 0) 

For larger x, the numerical studies£*2ii& generally find 
that both K and A s drop monotonically with increasing 
x. By the time x = X\ m 0.1, K is generally found to be 
close to 1, and by x — x c w 0.3, A s has dropped to values 
that are indistinguishable from 0, and K w 0.5. Thus, 
over most of the entire Luther-Emery phase, both the SC 
and the CDW susceptibilities are divergent. However, 
the SC susceptibility is the more divergent only at rather 
small values of x < x\. 

Although the charge fields exhaust the low energy de- 
grees of freedom of the Luther-Emery liquid, when we 
come to consider the effects of the single-particle hop- 
ping perturbation with small coupling constant St, we 
need to consider (as virtual intermediate states) high en- 
ergy states with the quantum numbers of an electron. 
Thus, we need to reintroduce gapped helds <p s to repre- 
sent the spin-degrees of freedom. Since this is standard 3 , 
we will not belabor the point; the appropriate continuum 
fermionic fields are 



±,a ~ exp • 



er#, ± ad J ± i 



Py/z} (7) 



where ± refer to left and right going fermions with mo- 
mentum near ±P/2, respectively, a = ±1 represents 
the spin polarization. It is important to stress that for 
strongly interacting problems, such as the present one, 



there is no simple relation between the original lattice 
fermions and the continuum fcrmion fields which describe 
the "physical " ^-fermions of Eq. 10. In particular, what 
appears as a 2kp CDW expressed in terms of ^-fermions, 
would be considered a4fcf CDW in terms of the original, 
lattice Fermions. In terms of these Vf-fields, the compo- 
nent of the charge density operator which varies with 
wave numbers near P is 

pp(v) = *I,<r*K,°- « ex P[* p y + iV2n<l>(y)] (8) 



while the singlet pair creation operator, 



1>(y) = [*£ * 



tx exp[i 



2tt6] 



(9) 



where in the right-most expressions we have again sup- 
pressed the dependence on the spin fields. 

Before leaving the single ladder problem, it is worth 
mentioning a useful intuitive caricature of its electronic 
properties. We picture a singlet pair of electrons on 
neighboring sites as being a hard-core bosonic "dimer." 
The undoped ladder can be thought of as a Mott insulat- 
ing state of these dimers, with one dimer per rung of the 
ladder, i. e. a "valence bond crystal" with lattice spacing 
one. To remove one electron from the system, we need to 
destroy one dimer and remove one electron, leaving be- 
hind a single electron with spin 1/2 and charge e. How- 
ever, when we remove a second electron from the system, 
we have the choice of either breaking another dimer, thus 
producing two quasiparticles with the quantum numbers 
of an electron, or of removing the unpaired electron left 
behind by the first removal, thus producing a new bo- 
son - a missing dimer - with charge 2e and spin 0. The 
persistence of the spin-gap upon doping the ladder can 
thus be interpreted as implying that the energy needed to 
break a dimer (of order A s ) is sufficiently large that one 
charge 2e boson costs less than two charge e quasiparti- 
cles. At finite x, the missing dimers can be treated as a 
dilute gas of hardcore bosons. That the elementary exci- 
tations of the undoped ladder can be constructed in this 
simple manner reflects the fact that this is a confining 
phas o 23 i 24 , not a spin liquid. 



III. INTER-LADDER INTERACTIONS 

We now address the effect of a small, but non-zero 
coupling (i e. single-particle hopping) between lad- 
ders, 5t > 0. Because of the spin-gap, St is an irrelevant 
perturbation in the renormalization group sense, and so 
does not directly affect the thermodynamic state of the 
system. However, second order processes result in var- 
ious induced interactions between neighboring ladders. 
These consist of marginal forward scattering interactions, 
which are negligible for small St, and potentially relevant 
Josephson tunnelling and back-scattering density-density 
interactions. 
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The important (possibly relevant) low energy pieces of 
these latter interactions are most naturally expressed in 
terms of the bosonic collective variables defined above: 

H'= -J2 [ dy{jcos[V^(6j - (10) 

+Vcos[(P j - P j+1 )y + V^ifa - 0, + i)]} , 

where Pj = 2nxj, with Xj the concentration of doped 
holes on ladder j, and 4>j and 8j are the charge field and 
its dual on each ladder. Here, again, the form of the low 
energy interactions between two Luther Emery liquids is 
entirely determined by symmetry considerations, but the 
magnitude of the Josephson coupling J and the induced 
interaction between CDW's, V, must be computed from 
microscopies; they are renormalized parameters which re- 
sult from "integrating" out the high energy degrees of 
freedom with energies between the bandwidth W ~ At 
and the renormalized cutoff, A s , or with wavelengths be- 
tween a and £ s = v s / A s where v s is the spin- wave veloc- 
ity. Thus, the dimensionless measure of the inter-ladder 
couplings, which for instance enter the expressions for T c , 
are J/W and V /W where W = A s /£ s . (So long as x is 
not too near x c , A s ~ J, and hence W ~ J.) 

Quantitative estimates of J and V could certainly be 
obtained, given the state of DMRG calculations, from 
studies of four-leg ladders consisting of two weakly cou- 
pled 2-leg ladders^. However, such calculations have 
not, yet, been carried out. Fortunately, our qualitative 
conclusions are not very sensitive to the values of V and 
3 , which can anyway be estimated with reasonable ac- 
curacy from bosonization, as discussed in Ref. [H The 
subtlety here is that the inter-ladder hopping is expressed 
in terms of microscopic lattice fermions, whereas our low 
energy theory is expressed in terms of the ^-fermions of 
Eq. J7J). However, since these have the same quantum 
numbers as an electron, and operate on the scale of A s , 
which is large with respect to St, there is no reason to 
expect any large renormalization of the hopping param- 
eters. If we assume that the inter-ladder hopping can be 
approximated as St times an operator representing the 
hopping amplitude for ^-fermions, then from second or- 
der perturbation theory we obtain 



J~V~A (St) 2 / J 
where A is the dimensionless function of A s /J 



(11) 



A = J Jdydr (e 



a i V V2[[e s (r)-9 B (O)+0 S (r)+0 3 (0)] » 



r = (y, t) and r denotes imaginary time, the expectation 
value is taken with respect to the spin-fields on the decou- 
pled ladders, and in deriving this expression we have as- 
sumed that the charge fields are slowly varying compared 
to the spin-fields. Simple scaling arguments of the sort 
discussed in Ref. suggest that 4 ~ 1 as A s /J ^ 0. 
(For further discussion see footnote 1261) In any case, so 



long as x is not too close to x c , A s is of order of the ex- 
change coupling, J, so it is reasonable to assume A ~ 1. 
The only aspects of this estimate which matter qualita- 
tively for our present purposes is that the two couplings 
are comparable in size, J ~ V and both are small in 
proportion to St 2 . 



IV. RENORMALIZATION-GROUP ANALYSIS 
AND INTER-LADDER MEAN FIELD THEORY 

The effect of these inter-chain couplings can be de- 
duced from an analysis of the lowest order perturbative 
renormalization group equations in powers of the cou- 
plings V and J . However, equivalent results are obtained 
from inter-ladder mean-field theoryi&il, which is concep- 
tually simpler. These equations are the analogue of the 
BCS gap equations applied to this model, and are ex- 
pected to give a quantitatively accurate estimate of T c for 
small St/A s for precisely the same reason. A discussion 
of the accuracy of interchain mean-field theory is given 
in Appendix [S] In the present two-dimensional system, 
T c should be interpreted as the onset of quasi-long range 
order, i. e. as a Kosterlitz-Thouless transition. 

To implement this mean-field theory, we need to com- 



pute the expectation value Mj(hj) = (cos[\j2irQj\) of the 
pair creation operator on an isolated ladder, where the 
expectation value is taken with respect to the mean-field 
Hamiltonian 



Hmf — Hj — hj 



J dy cos[\/2it9j] 



(12) 



in which Hj is the effective Hamiltonian in Eq. JBJ with 
parameters appropriate to ladder j, and hj represents the 
mean- field due to the neighboring ladders, and so satisfies 
the self-consistency condition, 



hj = J[M j+1 +Mj_ 1 }. 



(13) 



The expression for the mean-field transition temper- 
ature can be expressed in terms of the corresponding 

susceptibility, Xsc = 9Mj(h)/dh\h=o, which is related 
to the superconducting susceptibility in Eq. (JTJ by a 
proportionality constant which depends on the expecta- 
tion value of the spin-fields. In the case in which all the 
ladders are equivalent, this yields the implicit relation 
2Jxsc{T c ) — 1. For an alternating array of A and B 
type ladders, the expression for the superconducting T c 
is easily seen to be 



(2J) 2 xfMTc) xTc(Tc) 



~{B), 



1. 



(14) 



Notice that in the case in which the A and B type ladders 
are identical Eq. I|14f> reduces properly to the expression 
for equivalent ladders. The expression for \sc from Eq. 
(JTJ can be used to invert Eq. (|14|) to obtain the estimate 
for T c given in Eq. 
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The mean-field equations for the CDW order are ob- 
tained similarly. The expression for the transition tem- 
perature for CDW order with wave-vector P is 

(2V) 2 x ( cdw(P, T c) XcDwiP^c) =1 (15) 

where the notation is the obvious extension of that used 
in the superconducting case. The best ordering vector is 

that which maximizes T c . For P = Pa, Xcflw(^Ai^) 
diverges with decreasing temperature as in Eq. J2J, but 

Xcz)w(Pa,T) saturates to a finite, low temperature value 

when T~v e \P A -P B \. Thus, even if Xcdw( p a> T ) di " 
verges more strongly with decreasing temperature than 
Xsc> there are two divergent susceptibilities in the ex- 
pression for the superconducting T c , and only one for the 
CDW T c ; so long as the inequalities in Eq. J3J are satis- 
fied, the superconducting transition preempts the CDW 
transition! 



V. THE x -> LIMIT 

Since K — > 2 as x — ► 0, there is necessarily a regime 
of small x in which the superconducting susceptibility 
on the isolated ladder is more divergent than the CDW 
susceptibility. Here, in the presence of weak inter-ladder 
coupling, even the period 2 striped Hubbard model (i. e. 
with e = 0) is superconducting. However, care must be 
taken in this limit, since, as mentioned above, the range 
of energies over which H in Eq. JJjJ is applicable van- 
ishes in proportion to x 2 . Fortunately, a complementary 
treatment of the problem, which takes into account the 
additional terms, the ellipsis in Eq. ffijl. can be employed 
in this limit. The small x problem can be mapped onto a 
problem of dilute, hard-core charge 2e bosons (with con- 
centration x per rung) with an anisotropic dispersion, 
E(k) = tkl - Jcos[2k x ]. (The 2 reflects the ladder peri- 
odicity.) Consequently, for small x, 

T c w 2tt \plji xF(x) ~ \St\ x (16) 

where F(x) ~ 1/ lnln(l/x) is never far from 1, and the 
logarithm reflects^ the fact d = 2 is the marginal dimen- 
sion for Bose condensation. (This result is not substan- 
tially different for the period 4 striped Hubbard model, 
so long as e is not too large.) There is a complicated issue 
of order of limits when both St and x are small; roughly, 
we expect that T c will be determined by whichever ex- 
pression, Eq. I|14[) or Eq. . gives the higher T c , but 
with the understanding that xsc must be computed tak- 
ing into account the terms represented by the ellipsis in 
Eq. © which cause the susceptibility to vanish as x — > 0. 

The period 2 striped Hubbard or t— J model indeed has 
a superconducting phase at small x, because this phase 
is confined to rather small x < 0.1, where T c is small in 
proportion to both St and x. Moreover, this may still not 
be enough to establish a mechanism of high temperature 
superconductivity. The situation looks even worse when 



the effects of weak disorder — are considered - when the 
disorder strength is greater than the intra-ladder energy 
scale Ep — 2irtx 2 , it is unlikely that any sort of super- 
conducting coherence will survive. 

For an array of alternating ladders, the range of x 
for which superconductivity dominates is much extended. 
This means that the maximum T c is much greater, and 
the superconductivity much more robust to disorder— for 
the period 4 than the period 2 striped Hubbard model. 



VI. OPTIMAL DEGREE OF INHOMOGENEITY 
FOR SUPERCONDUCTIVITY 

Here we have established that in a strongly striped 
Hubbard model, superconductivity is produced directly 
by the repulsive interactions between electrons. The re- 
sulting T c is proportional to a positive power of St/t, and 
so rises as the stripe order becomes less strong. It is 
thus natural to ask: Is the stripe order introduced in the 
present paper simply a calculational crutch which per- 
mits us obtain well controlled results or is inhomogeneity 
essential to the mechanism of high temperature super- 
conductivity, as has been suggested&2£*22i22i2i in several 
previous studies? 

The answer to this question turns on the issue of 
whether or not the uniform Hubbard model, and its 
strong coupling relative the t — J-model, by themselves 
support high temperature superconductivity. This ques- 
tion has been the focus of much theoretical research since 
the discovery of superconductivity in the cuprates. To 
this date this is not a settled issue. Nor is it the purpose 
of the present paper to review this extensive literature. 
Variational calculations have been interpreted both as 
giving evidence in support^ and against superconduc- 
tivity in Hubbard and t — J-type models. There is also 
considerable evidence, from several numerical techniques 
and high temperature expansions, that the canonical t— J 
and Hubbard models on a square lattice most likely do 
not support high temperature superconductivity; instead 
they show clear evidence for other types of order which 
compete with superconductivity^^. 

Assuming that the uniform model does not support 
high-temperature superconductivity, it follows from the 
arguments given in the previous sections that there is an 
optimal degree of inhomogeneity (an optimal degree of 
stripe order) for a strongly correlated system to exhibit 
superconductivity. Probably this occurs when St ~ A s . 
An analogous result was established^ recently in the 
weakly interacting limit of the 4-leg ladder (itself a car- 
icature of a single unit cell of the present model). We 
should also note that there is nothing essential about hav- 
ing period 4. In fact, the longer the period the more the 
CDW instability is suppressed and the larger the range 
of superconductivity. 
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VII. RELATION TO SUPERCONDUCTIVITY 
IN THE CUPRATES 

While the main purpose of the present paper was to 
establish, as a point of principle, that the striped Hub- 
bard model analyzed here exhibits high temperature su- 
perconductivity, a few comments are in order concerning 
the more general implications of the present results for 
the mechanism of superconductivity in the cuprates. 

Firstly, the explicit striped inhomogeneities introduced 
here are a caricature of the spontaneous symmetry break- 
ing in a charge striped phase. However, the model pos- 
sesses a large spin-gap, and so does not contain any of the 
physics of low energy incommensurate spin-fluctuations 
which are the principle experimental signatures to date 
of stripe correlations in the cuprates. Secondly, although 
the superconducting state is "d-wave-like" in the sense 
that the order parameter changes sign under rotations 
by 7r/2, since the striped Hamiltonian explicitly breaks 
this symmetry, there is no precise symmetry distinction 
between d-wave and s-wave superconductivity. Thirdly, 
the superconducting state is not even truly adiabatically 
connected to the superconducting state observed in the 
cuprates, because the existence of a spin-gap implies the 
absence of gapless "nodal" quasiparticles in the super- 
conducting state. However, the transition between a 
node-less and nodal d-wave-like state was studied in Ref. 



35 36], where it was found to be a mean field (Lifshitz) 
transition with relatively little effect on T c , Moreover, 
using the same lines of reasoning employed in that arti- 
cle, it is possible to make compelling (although not en- 
tirely rigorous) arguments that upon heavier doping, the 
present model, too, will exhibit a nodal superconduct- 
ing state. We are currently working to obtain a more 
complete treatment of the phase diagram of the present 
model. 

The present model realizes the idea that the pairing 
scale, in this case the spin-gap, can be inherited from 
a parent Mott insulating state. Moreover, like the un- 
derdoped cuprates, the gap scale in the present model 
is a decreasing function of increasing x, while the actual 
superconducting transition occurs at a T c much smaller 
than A s /2, and is determined by the phase ordering tem- 
perature rather than the pairing scale. Hence, for x not 
too close to x c , this model exhibits a pseudogap regime 
for temperatures between T c and T* ~ A s /2, reminis- 
cent of that seen in underdoped cuprates. However, T c is 
always bounded from above by A s and so tends to zero 



as x 



Xr. 
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APPENDIX A: ACCURACY OF THE 
INTERCHAIN MEAN-FIELD THEORY 
ESTIMATES 

In this section, we discuss the accuracy of the inter- 
chain mean-field theory. Although no general proof ex- 
ists (to the best of our knowledge), we believe that it is 
asymptotically exact in the present case, at least to log- 
arithmic accuracy (as defined in Eq. (£H) below). The 
latter conclusion also follows, as mentioned in the text, 
by comparison with perturbative RG calculations. 

Quite generally, using an argument based on Grif- 
fiths inequalities, one knows that the exact T c of a gen- 
eral anisotropic ferromagnetic system (not necessarily 
an Ising model) will obey the bounds: T C (J X ) < T c < 
T c (J y ), for < J x < J y , where T C (J) is the T c of an 
isotropic system of coupling constant J. However, for 
specific systems it is possible to establish more precise 
estimates of T c . 

As our first example, consider the 2D anisotropic Ising 
model on a square lattice, with couplings J x and J y < J x 
in the x and y directions, respectively. In particular, for 
the case of the 2D Ising model, it is also known that the 
exact T c is the solution to the equation^ 



sinh(2 J x /T c ) sinh(2 J y /T c 



1. 



(Al) 



Interchain mean-field theory for the same model gives the 
familiar expression for the mean-field transition temper- 
ature, Tq: 



2J v Xid(T ) = 1 
which is analogous to Ea. (|14l) . and where 
Xid(T) =T- 1 cxp[2J x /T\ 



(A2) 
(A3) 



is the susceptibility of the ID Ising model. In the limit 
of small J y /J x , it thus follows that the ratio 



In 2 



T c \n[J x /J y 



+ 



(A4) 



tends to 1 as J y /J x — > 0, i.e. the interchain mean-field 
theory is asymptotically exact without any apologies. 

Before leaving the Ising example, it is interesting to 
see how well the interchain mean-field theory works when 
extrapolated to the isotropic case J x = J y = J. It is easy 
to verify that T = 3. 53 J and T c = 2 J/ ln[l + y/2], so 



T /T c = 1.55 for J y /J x = 1. 



(A5) 



In general, Tq/T c rises monotonically from 1 for increas- 
ing J v /J x , but To gives a reasonably good estimate of 
T c over the entire range of parameters. (Note, ordinary 
mean field theory gives Tq — 2[J X +J y ], which is not much 
worse than interchain mean-field theory in the isotropic 
limit, but Tq/T c — > oo as J y /J x — > 0.) 

Now, we move to the 2D classical XY model on a 
square lattice - the case of most direct relevance to the 
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estimates of T c made in the text. The susceptibility of 
an isolated chain can easily be seen to be 



Xid(T) = — 



1 fI a (J x /T)+h(J x /T) 



W x /T)-h{J x /T) 



(A6) 



where I n (x) is a Bessel function. For J y /J x <C 1, Eq. 
(|A3|) and Eq. (|A6fl yield the following estimate of the crit- 
ical temperature 



Tq — 2y J x J y 



l + O 



\J Jy/Ja 



(A7) 



while To = 1.755 J in the isotropic limit. 

Unlike the Ising case, no exact results exist for the 2D 
XY model. Extensive Monte-Carlo work has been done 
on the isotropic 2D XY model, form which we know2& 
that the Kosterlitz-Thouless transition occurs at T c = 
0.89J, so in this limit Tq/T c ss 2. In the limit of extreme 
anisotropy, the 2D classical XY model can be mapped 
onto the familiar ID quantum XY (rotor) model with 
Hamiltonian 



L' 2 



A 



— - - COS(tf„+i 



V) 



(A8) 



where the coupling constant is A = 2J x J y /T 2 (see Ref. 
[39|) The critical value of the coupling of the quantum 
rotor model, A c , has been computed quite accurately 
using a Padc-Borcl rcsummation of the strong-coupling 
series^S^i. Using the notation of these papers, an accu- 
rate estimate for the critical coupling to be A c = 1.8±0.5 
is obtained. By carefully inverting this mapping, we get 



T c — AyJ J X Jy 



l + O J Jy/ 7 : 



where A = 1.05 ± 0.1. Thus, we see that 
T Q /T C ^(2/A) as J y /J x - 



(A9) 



(A10) 



It seems unlikely that the error bars on A are sufficient to 
be consistent with a limit of 1 . The interchain mean-field 
theory is therefore found to be asymptotically exact only 
to logarithmic accuracy, i.e. 



In T 

In T c 



(All) 



This, we believe, is the generically true of interchain 
mean-field theory as applied in the present paper. None- 
the-less, in all cases where the exact answers are known, 
interchain mean-field theory gives estimates of T c that 
are within a factor of 2 of the exact results. This is cer- 
tainly sufficiently accurate for present purposes. 

Finally, it is worth mentioning that the 2D XY model 
is something of a worst-case example, because 2D is the 
lower critical dimension and hence fluctuation effects are 
anomalously large. If we consider an anisotropic 3D XY 
model with couplings J x > J y > J z in the three direc- 
tions, the mean-field transition temperature can still be 
readily computed according to 2(J y + J z )xid(To) = 1. 
Monte-Carlo results exisl^ for the T c of layered models, 
Jy = Jx = J for various values of J z /J x . For instance, 



T c = 1.1J T Q /T C = 1.60 for J z /J = 0.01 
T c = 1.324/ T /T c = 1.41 for J z /J = 0.1 
T c = 2.2 J T /T c = 1.29 for J z /J = 1.0. 



Clearly, even a very small amount of interplane coupling 
can be expected to greatly improve the accuracy of our 
T c estimates. (Interplane mean-field theory, of course, is 
still more accurate, as shown in Ref. |42^.1 
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